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Abstract 



The scaling limit of the critical percolation, is it a black noise? 
The answer depends on stability to perturbations concentrated along 
a line. This text, containing no proofs, reports experimental results 
that suggest the affirmative answer. 

Introduction 

By percolation I mean the two-dimensional critical site percolation on the 
triangular lattice (see for instance 3J). Its (full) scaling limit exists, is con- 
formally invariant and may be described by a random countable set of non- 
crossing SLE 6 curves (see 0). 

By a noise I mean a homogeneous continuous product of probability 
spaces, as defined by Def. 3dl]. The scaling limit of percolation could 
lead to a noise (and the noise should be black, as defined by Def. 7al]); 
this idea was discussed informally Question llbl] but encountered the 
following difficulty. 

Let a smooth domain D C R 2 be split by a smooth curve in two domains 
Di, Di- The random countable system of non-crossing SLEg curves in D, - 
call it the random configuration in D, — being restricted to D\ and D 2 gives 
two independent random configurations in D\,D2- Sewing them together 
appears to be a subtle matter. Does it admit any freedom? In other words, 
is the configuration in D uniquely determined by its restrictions to Di and 
D 2 1 

More formally: the limiting model on the whole plane R 2 is described 
by a probability space (fl,^, P) (of configurations). For every open interval 
(s, t) C R we consider the domain (s, t) x R c R 2 (strip) and the correspond- 
ing 'local' sub-a-field jF s t c T . Clearly, jF rs and !F Sjt are independent for 
r < s < t, and the upward continuity is satisfied: 

e>0 
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(that is, JF s t is the least sub-cr-field containing all J r s+£jt - e )- The question is, 
whether 

(1) rt t = ^F r ,s V T S} t 

or not; that is, whether the least sub-a-field T r jS V T s ^ containing jF r s and 
T s t t is the whole jF r t , or its proper sub-cx-field. The latter case would mean 
that the line {s} x M c I 2 cannot be neglected. Note that triviality of the 
cr-field J- s -, s + = n e> oJ? r s- e ,s+ £ is only necessary, since in general n e >o(.7>,s V 
Fss^+e V F 8} t) may exceed T r>a V (n £>0 .F g - e , s + £ ) V F s>t (see for detailed 
analysis). Note also that the difficulty cannot be avoided by waiving the 
upward continuity, since the latter holds for every noise Prop. 3d3]. We 
may reformulate (0) as 

(2) T Tlt = \/ {F r ,s-e V F s+e ,t) • 

£>0 

The question is, whether percolation (or rather, its scaling limit) is stable 
under a strong perturbation in an infinitesimal strip, or not. This strong 
concentrated perturbation is quite different from the distributed weak per- 
turbation examined in pQ. Sensitivity of percolation established there is 
micro-sensitivity in terms of [H Sect. 5c/5.3], where it is opposed to block 
sensitivity. The latter also corresponds to a distributed weak perturbation, 
but this time the correlation length of the perturbation tends to slower 
than the pitch of the lattice. Block sensitivity of percolation is noted in jH 
8a2/Remark 8.2]. 

Two possibilities remain open. Maybe, percolation is strip stable, that is, 
satisfies (0). Then it leads to a noise, and the noise is black due to the block 
sensitivity. Or maybe percolation is strip unstable (that is, violates 0). 
Then it does not lead to any noise. A trivial example of strip instability 
(and moreover, strip sensitivity) is given by the worst Boolean function, 
- the product of random signs (over all lattice points in the given domain). 
Here, trying scaling limit, we face an obstacle, — a continuum of independent 
random signs. 

1 The idea of the experiment 

Being unable to prove or disprove the equivalent properties (0), © for the 
scaling limit of percolation, I conducted an experiment as follows. First, 
the domain shown on Fig. 0a) (the union of two equilateral triangles) is 
filled in with percolation data. That is, each hexagon gets black or white 
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with probabilities 0.5, 0.5 independently of others. However, the boundary 
is colored deterministically, its left part in white, right part in black. The 
corresponding exploration path is constructed. Second, the percolation data 
near the 'equator' are removed and generated anew. The second exploration 
path is constructed. Finally, the two exploration paths are compared. This 
is a single trial. Many such trials are conducted, independently of each other, 
for different values of parameters (see Sect.Ej). 




(a) (b) (c) (d) 

Figure 1: the domain (a); percolation data and the first exploration path (b); the 
percolation data near the 'equator' are removed (c) and generated anew, giving 
rise to the second exploration path (d). 

The comparison of the two exploration paths 71,72 is based on the fol- 
lowing metric borrowed from [21 eq. (2)]: 

(3) ^(71,72) = mfmffi |ti (0 - 72(01 1 

here \ ji(t) —72(01 is the usual Euclidean distance between the two points of 
the plane (the unit of length is such that the diameter of our domain equals 
2, that is, each side of the rhombus equals 2/y / 3), and the infimum is taken 
over all parameterizations of the paths. Similarly to , a path is treated as 
an equivalence class of continuous functions [0, 1] — > M. 2 modulo monotonic 
(increasing) re-parameterizations. 

2 Experimental results 

Each trial is described by two parameters n, k as follows. The number of 
hexagons on the 'equator' (including two boundary hexagons) equals n + 2; 



3 



that is, the (discrete) rhombus consists of 2n + 1 rows (of 2, 3, . . . , n + 1, 
n + 2, n + 1, . . . , 3, 2 hexagons respectively). Out of these In + 1 rows, k 
rows are removed and generated anew. The set of k rows is either symmetric 
w.r.t. the 'equator' (if k is odd), or contains one additional row above the 
'equator' (if k is even). For instance, n = 8 and k = 2 on Fig. ^b,c,d). 

Table ^ represents 27 samples, for 27 pairs (n, k), of 250 trials each. For 
every sample, the median of 250 values of the distance is reported. 



^\ n 
k 


16 32 


64 


128 


256 


512 


1024 


1 


0.22 0.19 


0.16 


0.13 


0.11 


0.09 


0.08 


2 


0.24 


0.18 


0.19 


0.15 


0.11 


0.10 


4 




0.27 


0.24 


0.19 


0.15 


0.12 


8 






0.28 


0.23 


0.19 


0.16 


16 








0.29 


0.22 


0.18 


32 










0.31 


0.22 



Table 1: Median values of the distance between the two exploration paths. 

The distance is basically defined by (j3J). However, for speeding up the 
computation, I use an approximation to (described in Sect. EJ), which 
leads to a systematic error bounded by ±0.03. Thus, each median in Table 
Qhas a systematic error bounded by ±0.03, and a sampling error, the mean 
square deviation being about 7% of the median. 

The results suggest that the typical distance tends to as the strip width 
e = - tends to 0. Maybe we observe a power law e a for some exponent a 
not far from 1/3. 



3 Examples 

The median (= 0.22) of the first sample (n = 16, k = 1) resulted from a 
definite trial (no. 21, in fact), — call it the median trial, — presented on 
Fig. El Only the first exploration path is shown on Fig. Efa) together with 
the hexagons relevant to it. On Fig. [2Kb) the path is not shown explicitly, but 
you can restore it easily. On Fig. Etc) the second exploration path is shown 
together with the hexagons relevant to the first path. Hexagons relevant to 
the second path are not shown explicitly, but you can restore them easily. 
Fig. Efd) presents the most interesting, 'equatorial' part of Fig. Etc). 
More examples follow (Fig. El EJ), in the same format as Fig. Efd). 
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(a) (b) (c) 




(d) 

Figure 2: The median trial for n = 16, k = 1. The first exploration path (a) 
and the hexagons relevant to it (a,b); the same hexagons together with the second 
exploration path (c,d). 




k = 1 k = 2 

Figure 3: Median trials for n = 32. 
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k = 4 

Figure 4: Median trials for n = 64. 
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4 Technicalities 



First, about the metric (jHJ) and its discrete approximation used. Treating the 
exploration path as a finite sequence of points, — vertices of the polygonal 
line (rather than the line itself), I define a distance d(a, b) between two such 
sequences a = (a , . . . , clm) and b = (b , . . . , b N ) as follows: 

d(a, b) = min max |a Mfc — b Vk \ ; 

(ji,u) k=0,...,M+N 

here \a^ k — b Uk \ is the Euclidean distance between the two points, and (/i, v) 
runs over all pairs of sequences /i , • • • , ^m+n and u , . . . , i>m+n of integers 
satisfying 

= no < - ■ ■ < Hm+n = M , = v < ■ - ■ < i/m+n = N , 
[ik + v k = k for k = 0, . . . , M + N . 

Clearly it is an approximation to (jHJ). (Not a metric, since d(a, a) ^ 0, which 
is harmless.) Dynamic programming helps us to compute the minimum over 
(/i, v\ To this end, introduce 

dk ti = min max |a w — b n \ 

{li,v);n k =i l=0,...,k 

for i = 0, . . . , min(fc, M) and consider the array D k = (dk,i)i=o,...,min(k,M)- 
Having D k we can compute D k+ i easily. Thus, we compute D , then D x , 
and so on, up to D M+N ; the latter gives us d(a, b). 

However, M + N is typically as large as 10 5 (for n = 1024, my worst 
case). For speeding up the computation, I replace the long sequence a = 
(ao, ... , om) with a subsequence (a kl , . . . , a kp ) such that \a ki — ai\ < e for 
ki < I < k i+ i. (The same is done for the other sequence b.) The parameter 
e was chosen to be 0.03. 

Now, other technicalities. Random colors (of the hexagons) are produced 
by a Wichmann-Hill generator of pseudo-random numbers, initialized (before 
each trial) with a one-to-one function of log 2 n, log 2 k and the number of the 
trial (within the sample). 

The program, in the 'Python' programming language, is available on my 
Web site. 

The computer: a personal computer with matherboard Intel D815EGEW, 
processor Intel Pentium III (933 MHz), cache RAM 256 Kb and SDRAM 256 
Mb (133 MHz). 

The operating system: SuSE Linux 9.0. 

Python: version 2.3+. 

One trial takes typically from 0.5 min to 1 min for n = 1024, and 10-15 
sec for n = 512. 
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Conclusion 



The results (Table QJ) suggest that the scaling limit of percolation is strip 
stable, and therefore may be treated as a black noise. 

One could object that the strong concentrated perturbation is not enough; 
a distributed 'microscopic' perturbation should be added, such that explo- 
ration paths remain macroscopically close to the same SLE 6 curves, but 
are changed microscopically. That is correct, but I guess that such 'micro' 
perturbations may be absorbed somehow by the strong perturbation in the 
infinitesimal strip. 
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